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ABSTRACT 



Context. Stellar and galactic systems are objects in dynamical equilibrium that are composed of 
ordinary baryonic matter hypothetically embedded in extended dominant dark matter halos. 
Aims. Our aim is to investigate the scaling relations and dissipational features of these objects 
over a wide range of their properties, taking the dynamical influence of the dark matter compo- 
nent into account. 

Methods. We study the physical properties of these self-gravitating systems using the two- 
component virial theorem in conjunction with data that embrace a wide range of astrophysical 
systems. 

Results. We find that the scaling relations defined by the properties of these objects admit a dark- 
to-luminous density ratio parameter as a natural requirement in this framework. We also probe 
dissipational effects on the fundamental surface defined by the two-component virial theorem and 
discuss their relations with respect to the region devoid of objects in the data distribution. 
Conclusions. Our results indicate complementary contributions of dissipation and dark matter to 
the orign of scaling relations in astrophysical systems. 

Key words, scaling relations - fundamental plane - elliptical galaxies - galaxy clusters - virial- 
ized systems - dark matter halos 



1. Introduction 

The present cosmological paradigm of structure formation relies on indirect evidence that a dark 
matter component dominates the dynamics of large structures (e.g. Padmanabhan 2006 , and refer- 
ences therein). Gravitational instabilities increase the initially small amplitude fluctuations in the 
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primordial universe, resulting in the evolution towards virialized systems like galaxies and clusters 
of galaxies (Peebles 1980 Coles & Lucchin|[l995 Dodelson 2003 ). However, this rather simple 
mechanism is in reality exceedingly more complicated, mainly because of non-linear effects, such 
as the role of dissipation of the baryonic component inside dark matter halos (Silk & Norman 1981 ) 
and the nature of dark matter itself (e.g. |Gao & Theuns|2007| l. In addition, in the hierarchical sce- 
nario of structure formation, low-mass objects are formed first, and then larger systems are formed 
by subsequent merging of smaller subunits, resulting to a greater or lesser degree in a dynamically 



complicated history of fusions for probably any given selfgravitating object seen today (e.g., De 
|Lucia & Blaizot|2007l|McIntosh et al.|2008) . 

Despite the evident difficulties reaching complete understanding of these mechanisms, it is 
remarkable that virialized systems do present some kind of global regularity in their scaling re- 
lations. For instance, it is found that all virialized stellar systems are reasonably confined to a 
two-dimensional sheet in their parameter space defined by three independent variables, derivable 



from quantities such as luminosity, radius, and projected velocity dispersion (e.g. Schaeffer et al. 



1 993 1 [Burstein et al.|1997| l. The study of elliptical galaxies, for instance, is of paramount impor- 
tance in this respect, given that their 2-D manifold of physical observables are a very well-defined 
plane (the so-called 'fundamental plane", c.f. Djorgovski & Davis 1987; Dress ier et al.|1987). with 
a characteristic small scatter and significant deviation ("tilt") from the simple virial expectation 
(see Tortora et al. 2009 and references therein). 

The idea of the fundamental plane (FP) has its origin in some restricted relations between global 
observables (central velocity dispersion, effective radius, and mean effective surface brightness) 
in stellar systems that emerge after statistical analysis of extensive photometric and spectral data 
accumulated for a number of objects over the past two decades (Djorgovski & Davis |1987||Dressler] 
et al.|1987||Djorgovski & de Carvamo[T990]|Bernardi et al.|2003[ ). The FP is a bivariate family of 



correlations that are well established for elliptical galaxies in the optical and infrared ranges (e.g. 



Guzman et al. 1993 Pahre et al. 1995 respectively). Other FP-like relations are observed for spiral 



and dwarf galaxies (e.g. Jabl onka et al.|1996}|Cody et al.|2009) > and galaxy clusters (Schaeffer et al.| 



1993 Fritsch & Buchert 1999 Lanzoni et al. 2004 Araya-Melo et al. 2009 1. Indeed, Burstein et al. 



( 1997| l propose the term "cosmic metaplane" to describe this ensemble of FPs from globular to 
galactic cluster scale. 

The existence of such a metaplane reinforces the idea of some dynamical regularity in the 
properties of a wide range of different astrophysical objects, which certainly deserves deep analysis 
of its own. But the cosmic metaplane also reveals some other interesting features of the structure 
formation process, such as a region that seems to be completely devoid of objects, called the "zone 
of exclusion" (ZOE, e.g. Bur stein et al.| 1 997 ). Such a feature should be accounted for by the current 
paradigm of galaxy formation and clustering in the Universe. 

|Dantas et al.| ( [20b 0) (hereafter Paper I) shed some light on these questions by putting forward a 
model based on the hypothesis that self-gravitating stellar systems have extended dark matter halos 
and described their equilibrium state by a modified, two-component virial theorem (2VT). They 
have found a fundamental surface in the parameter space that embraces the properties of a wide 
range of stellar and galactic systems. This surface shows remarkable compatibility with the cosmic 
metaplane discussed by Burstein et al. ( 1997\ and could partially explain its existence from simple 
dynamical arguments. 
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It turns out that the 2VT hypothesis, however compelling it may seem considering the wide 
scope of compatiblity resulting from few fundamental assumptions, does not necessarily constrain 
or contradict the few other alternative explanations to the FP considered in the literature (see, 
e.g., Proctor et al.|2008 Tortora et al.|2 009 and references therein). Indeed, the nature of scaling 



relations such as the FP is still controversial at this time. In order to advance our understanding 
of the 2VT hypothesis as the main factor in explaining the scaling relations of virialized systems, 
we analyze possible paths to disentangling additional effects from the purely dynamical one as 
imprinted in the 2VT relations themselves. 

A usual and important additional nondynamical effect to explain the FP is dissipation. For 
instance, Robe rtson et al.| ( p006[ ) demonstrate that if disk galaxy mergers are strongly dissipational, 
tilt is generated in an FP, in agreement with observational determinations. Actually, some authors 
consider dissipation as a necessary and sufficient condition for explaining the nature of the FP 
(e.g. Hopkins et al. 2008 ). This possibility motivates an analysis of how far the 2VT hypothesis is 



orthogonal to the dissipational hypothesis, if at all. For that we will define quantitative estimators 
for dissipation. The most interesting features to analyze under these estimators is the ZOE, a subject 
that has not been addressed much in the past ~ 10 years since the publication of Paper I (see Zaritsky 
et al.|2006) Gadotti|2009 for recent examples of the ZOE for spheroids and bulges, respectively). 



However, in that work, the origin of the ZOE and the gap between stellar and galactic systems was 
not addressed, so this analysis is offered here. 

This paper is organized as follows. In the next section, we review the 2VT hypothesis and 
possible evidence for it. We outline the interplay between dissipation, the 2VT, and the ZOE in 
Section[3] Finally, we discuss our results (Section|4]i. 

2. Dark matter: scaling relations and 2VT 

In this section we revisit the 2VT formulation and study the evidence for it based on a statistical 
analysis of the 2VT fittings to a wide range of astrophysical objects. 

2.1. Revisiting the 2VT 

The properties of regular nearby galaxies are slowly evolving at present, after a presumably pre- 
vious and complex gas-rich protogalactic stage. Thus, observables such as surface brightness, ef- 
fective radius, mean velocity dispersion, and luminosity can be used to specify the parameters of 
these objects. At the same time, when considering galaxies as stationary self-gravitating systems, 
their physical properties should at some level reflect an equilibrium state through the virial relation, 
r e oc o-q/^ 1 (using the one-component virial theorem, 1VT), which implies that only two of a great 
number of observables are independent. Actually, the FP describes a global relation of virial type, 
r e oc cr£/f but, differently from 1VT, with exponents A ~ 1.53 and B 0.79 for elliptical galax- 
ies observed in the near-infrared (Pahr e et al.p 998). This means that there is no direct mapping 



between the virial plane and the FP. The reason for the difference of these relations (also called the 
FP tilt) is a question of current debate involving different hypotheses such as the non-homology of 
ellipticals properties or a mass-luminosity ratio variable with the total luminosity (see, for instance, 
Capelato et aL]T995l \TWf\ |Hjorth & Madsen|1995||Ciotti et al.|1996[|Kritsuk|1997l|Bekki|1998| >. 



Trujillo et al. (2004) suggests that the FP tilt could be a combination of stellar population plus 
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non-homology effects, while [Lanzoni et aT7| ( |2005) > points out scaling relations of ellipticals might 
be the combined result of cosmological collapse plus dissipative merging. 

When extending the problem to other scales and systems, we find a similar situation: the cos- 
mic virial plane described by the 1 VT, defined as the ensemble of all possible collapsed objects of 
all masses and radii (Burstei n et al.|1997) , does not coincide with the observationally determined 
cosmic metaplane, as shown in Paper I. This metaplane is also tilted with respect to the virial 
expectations. This problem is mitigated in Paper I by adding a term to the virial equilibrium equa- 
tion considering the gravitational energy due to the interaction of dark and baryonic components. 
Applying the 2VT to stellar and galaxy systems, they find 

2 GM LUM An 2 

{v > = + —Gp DM {r ) LUM , (1) 

rwM j 

where pom is the mean density of the dark matter halo within the region containing the luminous 
component, (v 2 )lum is the mean square velocity of the stars (or galaxies) and 

2 _ fr 2 p W M(r)dV 

(r )lum = —j. , (2) 

J PLUMVndV 

where plum is the baryonic density inside tujm- In terms of observed quantities, Eq.[T]amounts to 

cr 2 = C*(I e r e + br 2 ), (3) 

in which <tq is the central velocity dispersion of the baryonic component, C* = 2nGC r C v (x) L 
and b = ~ (t) lum Pdm> with {t) lum the bar y° nic mass-to-light ratio, and R = {r 2 ) LUM /r 2 . For 
Pdm = 0, we recover the 1VT, cr 2 = 2nC r C v GM e / 'r e , where M e = M tot /2. 

The kinematical-structural coefficients (C r , C v ) translates the physical into observational quan- 
tities, and they may or may not be constant among galaxies: cr^ = C,.(v 2 ) (where (v 2 ) is the mean 
square velocity of the particles), r e = C r r c (where r c is the gravitational radius of the system). In 
this work, we assume R/C r ~ 20 and C r C v * 0.2 as typical values for these parameters (see Paper 
I). 



2.2. Replotting the 2VT 



In Paper I, we worked in the k parameter space (see Bender et al.|1992 1, in which the 2VT can be 
expressed as 



K3 



logC* 



1 



V3 V3 
where 

a) = (Ki - V3k 2 )/ V2 



log(l +M(T), 



log 



(4) 



(5) 



and thus u> describes the central luminosity density of the stellar systems. In this space, the tilt can 
be understood as a systematic increase in k^, along the FP (e.g. Ciotti et al.|1996| l. Also, from Eq. 
[4j we see that depends on C, b, and a>. 

Now, working with directly observable quantities and recalling that the total mass-to-light ratio 
(refered to B-band luminosities) within the effective radius r e is given by log [M e /L e ](Q) = V3/C3 - 
log C e ff (see Burstein et al.|1997 l, we can rewrite the 2VT as 



log 



= log C* + log (1 + MO") - log C eff 



(6) 
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or simply 



M, 



C" 



C, 



eff 



1 + 



PL 



(7) 



where C e /f = 2nGC r C v and we have defined pi = l e jr e and sow = logp L . 

For pom = 0, the only way of increasing (M e /L e ) along the FP is through a systematic variation 



in C r and C v , which would imply that galaxies form a nonhomologous family of objects (Graham 



& Colless|1997 Capelato et al.|1993 1997 1. Another possibility for increasing C* is to suppose a 
systematic trend in the mass-to-light ratio with galaxy mass: (M/L)lum M a (e.g. Dressier et al. 



1987| >. However, if p^M + 0, we can retrieve homology and still increase M e /L e (or K3), depending 



on the variation in the ratio b/pi, which is the dark-to-luminous density ratio parameter of the 2VT 
formulation. 




log a (km s"') 

Fig. 1. Velocity dispersion versus 0.5 log [C*(I e r e + br^)] (in km s _1 ). The solid lines represent 
the 2VT constrained by C* = 8.28 and p DM = 2.3 x 1(T 2 M G pc~ 3 (b = 0.2) (stellar systems) or 
Pdm = 5.8xl0~ 6 M Q pc~ 3 ) (b = 4.8xl0~ 5 ) (galactic systems). The symbols are for globular clusters 
(open squares), spiral galaxies (filled squares), elliptical galaxies (filled circles), groups dominated 
by elliptical galaxies (open triangles) and clusters of galaxies (filled triangles). The 1 VT prediction 
is equivalent to the 2VT line for stellar systems. 



In Paper I, a brief analysis was given for the adjusting procedure, which relies on certain as- 
sumptions about the mass-luminosity relation of the baryonic component, as well as on considera- 
tions about the DM density profile and spatial scalings for the baryonic and non-baryonic compo- 
nents. We refer the reader to that paper for the details. In particular, the values of pdm are found after 
adjusting for the parameter b, which gives two "typical" values for a good global fit to the dataset, 
according to the type of system: globular clusters and galaxies result in pom ~ 2.3 x 10~ 2 M pc -3 , 
and galaxy groups and clusters favor pom ~ 5.8 x 10~ 6 M Q pc~ 3 . For the present work, we as- 
sume the same adjusting assumptions (and explore the rationale behind this in a moment). Also 
in accordance with Paper I, we keep C e /f = 6.76 and C* = 8.28 constant (which corresponds to 
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Fig. 2. Effective mass-to-light ratio (in M G /L Q units) as a function of luminosity density (L Q pc -3 ). 
The solid lines represent the 2VT constrained by C* = 8.28, C e tf = 6.76, and p DM = 2.3 X 
1(T 2 M pc~ 3 (b = 0.2) (globular clusters and galaxies) orp DM = 5.8xl(T 6 M Q pc~ 3 (b = 4.8xl0~ 5 ) 
(galaxy groups and clusters), while the slashed line is the 1 VT. 

(M/L)lum ~ 1-6) for the data fittings. Based on these parameters, we plot in Figs. [T] and [2] the 
scaling relations defined by Eqs. [3] and [7] with the observational data provided by |Burstein et al. 
1997| l. We clearly see that, in the region corresponding to high central density luminosities, both 



1 VT and 2VT fit the points. However, in Fig. [T] as <x increases, only the 2VT relation is able to 
fit galaxy systems. At the same time, In Fig. [2] as pi decreases, only the 2VT relation fits all the 
observational data. 

In the next two sections, we analyze whether the assumptions on pom (namely, the assumed 
bi-modality in the dataset) are valid at the level of investigation that we aim at in the present work. 



We first (Section 2.3 1 quantify the evidence for the 2VT model vis-a-vis with the 1 VT model, with 



the given adjustment assumptions above, using statistical arguments. Subsequently (Section 2.4 1 



we perform statistical pairwise comparisons in order to group systems by relative affinity in their 
scaling relations. 

2.3. Evidence for 2VT 



To compare 1VT and 2VT relations, we used the Akaike information criterion (Akaike 1974). 



The differences in AIC values needed for model comparison only depended on the ratio of the 
maximized model likelihoods and the difference in the number of parameters between the two 
models (e.g. |Burnham & Anderson 2002): 



A/Ci - AIC 2 



-21a|^l + 2(&i-& 2 ). 



(8) 



Royall ( 1997 1 suggests that likelihood ratios equal to 8 and 32 are equivalent to setting significance 
levels of 5% and 1%, or AIC differences of 2 and 3.5, respectively. That is, when we have AAIC > 
2, the two models are significantly different. In the present analysis, we found that data reveal a 
compelling difference bewteen the models, with strong evidence of 2VT. Globular clusters and 
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galaxies have A(AIC) = 2.35 and A(AIC) = 45.90 for the scaling relations ^ and 0. At the 
same time, galaxy groups and clusters have A(AIC) = 34.31 and A(AIC) = 18.78 for Q and (|7]), 
respectively. 



2.4. Multiple comparisons 

It is important to ascertain whether, at the present level, a refinement in the adjustments of the 
2VT to the several types of systems analyzed here is an essential prerequisite for the subsequent 
formalism (dissipation parameter in the 2VT) that we introduce. 

Indeed, it would be desirable to have a "minimalist" model to explain the "cosmic metaplane" 
as a whole, with the least possible number of free parameters, and this was in fact the initial mo- 
tivation for the 2VT formulation. It became clear, however (see also Figs. [T] and [2}, that at least 
two "typical" values of p^M were needed under simple assumptions. Evidently, a more thorough 
investigation would have to address how different assumptions on the already mentioned mass- 
luminosity relation of the baryonic component, the DM density profile, and spatial scalings for the 
baryonic and non-baryonic components affect the 2VT fittings, hence the pr,M values. 

The most important issue would be the strength of the main 2VT result at present, namely, that 
the Pom f° r galaxies is about 3 orders of magnitude greater than that of groups and clusters of 
galaxies. It is not our intention in the present work to study the various fine-detailed adjustment 
values of p DM for these systems, although we acknowledge the expectation that their correspond- 
ing values could vary systematically with mass (e.g., Bullock et al.|200T "p] What we aim at here 
is simply given fixed baryonic and non-baryonic density profiles, show that it is possible to formu- 
late a further unique quantity in terms of the 2VT variables that encode dissipational information. 
Therefore, for the sake of our argumentation, we fix our adjustment assumptions to match those of 
Paper I, but acknowledge here that scaling variables, such as (r 2 )/,f/M, are actually functionals of 
the density profile assumed. 

To certify how possible conclusions derived from our formalism could be significantly affected 
by the fixed assumptions, we must analyze whether the bi-modality of pum makes sense statisti- 
cally. To answer this question, we present a pairwise comparison of means related to the variables 
included in relations ([3]l and Q. We used the Tukey "honest significant difference" technique 
(Tukey's HSD), a quite conservative approach to getting information on the pattern of differences 
between the means of specific groups of objects (e.g. Hayter|1986| l. We aim here to assure whether 
it is statistically significant to group by affinity galaxies, on one hand, and groups/clusters, on the 
other, from the given scaling relations of these systems, without reference to their possible specific 
Pdm values. 

Tukey's HSD test is a method of ensuring that the chance of finding a significant difference 
in any comparison is maintained at the confidence level of the test. For the present test, therefore, 
pairs of objects related by affinity in their scaling relations will have their comparison values com- 
patible with zero. In Fig. [3] we see the 95% confidence-level Tukey pairwise comparisons from 
our data. In upper and lower panels we plot the differences in the variables involved in expressions 
([3]) and Q, respectively. The confidence intervals reveal that some differences are not significant 
among two or more families of objects per variable. These internal similarities probably set the two 



We aim to investigate that issue in a future work. 
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distinct families that lead us to the tracks in Figs. [TJ and [2] under the assumption of two "typical" 
Pom values. In other words, pairwise comparisons allow us to conclude that the double-adjustment 
found in Section [272] is statistically acceptable. Therefore, for our present purposes, we regard our 
analysis as valid up to the point that (i) any scaling variables in the formalism are functions of the 
density profile, and (ii) fine-graining adjustments of the 2VT for individual systems under fixed 
assumptions of the density profile are not fundamental for the formalism. 




Fig. 3. Tukey pairwise comparisons for the scaling relations (|3j and (7j) at 95% confidence 
level. Abbreviations for astrophysical systems: GC (globular clusters), S (spiral galaxies), E (el- 
liptical galaxies), Gr (galaxy groups), and CI (galaxy clusters). In the upper left panel, we use 
C* = 2nGC r C v 1^) , instead of a constant value. 

' V \L )LUM 



3. Dissipation formalism, 2VT, and ZOE 

While galaxies are examples of baryonic cooling and strong dissipation, groups and clusters have 
too long a cooling time to have dissipated a significant amount of energy on their own scale (e.g. 
Blum enthal et al.|1984 1. The virial theorem (both 1VT and 2VT) cannot directly access the details 



of the dissipation mechanisms during structure formation. However, a slight reformulation of the 
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2VT allows us to build a new quantity to measure the luminous mass surface density of the self- 
gravitating systems, Elum> by relating two independent FP parameters, <xq and r e : 



1 o"o ^Pdm (r 2 )LUM 



(9) 

On dimensional grounds, this quantity may be relevant in considerations of the role of dissipation 
in these systems. 

In Fig.|4j we plot Elc/m versus ctq using data from Burstein et al. ( 1997 1. This projection should 
be visualized keeping in mind that it actually represents a projection of a higher (3D) dimensional 
parameter space (with r e being the third axis). Straight lines indicate the fits to the 1VT relation, 
whereas 2VT fits are characterized by curves that separate at some distinct points from the 1VT 
lines. The characteristic scales r e (assumed to be the median of the objects in each kind of system) 
are indicated as labels on the diagrams. 
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Fig. 4. Luminous mass surface density against velocity dispersion. The solid and dashed lines are 
the 2VT and 1 VT, respectively, for different values of r e (kpc). 



The most distinct feature of this figure is that, except for globular clusters, the 1VT implies 
a bad fit to data, unless one is willing to adjust an arbitrary series of 1VT relations for different 
r e scales in order to cover all data. On the other hand, the 2VT curves naturally bend along the 
2n/M-axis and follow the data distribution, making therefore a more "economical" way to describe 
these systems. In particular, for the case of galaxy groups and clusters, typical values of r e require 
the 2VT bend to fit the points, while the 1 VT lines would imply too large a r e for these systems. 

We should also point out in this figure that there are regions defined by lines that are almost 
perpendicular to the 1VT lines above, and below them there are no objects at all. For instance, the 



so-called "zone of exclusion" (ZOE) pointed out by Burstein et al. ( 1997 1 indicates an absence of 
high [cro, 'Llum] objects in the Universe. Correspondingly, there is another limit to these quantities 
precluding the existence of very low [cro, 'Llum] objects. The 2VT somewhat constrains this behav- 
ior since the 2VT curves tend to go asymptotically towards higher cr for a given r e (in comparison 
to the 1VT case). A remaining question is to know how to use the virial equations to estimate the 
depletion of objects in certain regions of the diagram, in particular the high [ctq, ~Llum\ corner. 
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According to Burst ein et al.| (fl997), the equation to describe the region allowed for objects is 

j e < 2.95 x 10 14 M; 4/ \ (10) 

where the effective luminosity density is j e = 0.75 x 10~ 3 I e /r e L pc~ 3 and logM e = log (cr^) + 
5.67. Recalling that a> = - log (I e /r e ), we have 

log (S^cr 4 ) < 0.75w + 7.52 + log G. (11) 

We define a new quantity to estimate the relative importance of dissipation from the equilibrium 
configuration of already formed astrophysical objects as 



Q = 







1 b) 




[Ceff) 


[pL p z L ) 



(12) 



The idea is that, as dissipation plays less of a role, p L decreases while M e /L e should increase, such 
that the larger Q, the less dissipation that occurs; therefore, the variation in Q with pi depicts to 
what degree dissipation dominates the formation process. Note that Q also depends on the dark- 
to-luminous density ratio parameter b/pi. Thus, using Q and recalling that E e = M e /^r^ and 
L e = nP e p L , we can use use expressions (7) and (1 1) to find 

< 

C 2 eff PlQri 

At the same time, taking a> = logp L ' and introducing (13) into (1 1), we get 
1 

Crfpi 



with C = C^~/(3.16 x 10 7 G). From this we can define 

, \ 0.36 



\ QminC r e I 

as the approximate maximum luminosity density for a given scale r e . 



The Q m i n values for the Burstein et al. (1997) sample are Q s,ar ~ 10 (globular clusters and 



galaxies) and Q*. ~ 10 5 M Lg 2 pc 3 (groups and clusters). Using these values in (19), we plot 
both p™" and the data points on the [r e ,pi\ plane (see Fig. [5]). In this figure, the Q lines indicate 
the behavior of the maximum dissipation efficiency (p™ flv ) allowed for a given scale r e . Globular 
clusters and galaxies have a linear relation between the two lines, while galaxy systems present a 
faster decrease of luminosity density for larger radii. In fact, two values of Q mm are necessary to 
regulate the distribution of points and prevent clusters and groups from settling in the ZOE region. 
This result is consistent with the idea of dissipation as the main vector organizing astrophysical ob- 
jects. However, dissipation and dark matter are tightly connected (see Eq. (7)), and we necessarily 
have different values of central dark matter halo densities when we pass from galaxies to clusters 
scales. In fact, if one ignores the dark component, dissipation can be significantly overestimated at 
cluster scales; hence, our result suggests some influence of pdm on the dissipative infall of baryons. 
At the same time, the infall itself can affect the dark matter component as well. As suggested by 
Blu menthal et al.| (fl986), dissipative infall can produce smaller dark core radii and higher dark 
central densities. This points out the idea of a deep connection between dissipation and the final 
configuration of both luminous and dark distributions of matter. 
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Fig. 5. Luminosity density as a function of r e . The solid lines indicate p™ ax for two values of Q, 
Dotted lines mark different locii in the plane. 



4. Discussion 

The virial metaplane with all its features encloses important aspects of galaxy formation. The 
phenomenology involved is rich and complex. In this work, we have concentrated our investigation 
on the 2VT formulation as a function of direct observables and the dissipational features on the 
fundamental surface defined by the 2VT as a continuation of Paper I. Our main results are 

- Working with directly observable quantities, we fitted the 2VT scaling relations, presented in 
Figs.[T]and[2] for a wide range of astrophysical objects. Statistical evidence based on AIC dif- 
ferences and Tukey's HSD test allowed us to conclude in favor of the bi-adjust 2VT description. 
This result indicates the importance of dark matter to the final configuration of astrophysical 
objects. 

- Using the dissipational parameter Llum, we found two regions devoid of objects, corresponding 
to the upper right and lower left corners on the [ctq; Y*lum\ plane (Fig. [4}. The non-existence of 
low [(To; "Elum] systems is probably related to clouds with too low temperatures and densities 
to form luminous objects (e.g. Schaeffer & Silk 1985), whereas the underpopulation of high 
[erg; I.lum] objects seems to agree with the work of Silk & Norman ( 1981| > in the context of gas 



clouds collisions, which indicates a limit for the kinetic energy being lost via dissipational clus- 
tering. The results for "Zlum (ZOE) suggest that dissipation probably plays the most important 
role in defining such features on the fundamental surface. 
- Defining a new dissipational estimator, Q, we found an expression for the maximum luminosity 
density as a function of r e and Q m j„ for each class of objects (c.f. Eq. 15). Once more, the ZOE 
is a clear feature in the [r e \pi\ plane. Since b can be written as a function of pi and M e /L e , this 
suggests a deep connection between dissipation and the final configuration of both luminous 
and dark distributions of matter. 



Taken together, these results indicate complementary contributions of dissipation and dark mat- 
ter to the orign of scaling relations in astrophysical systems. This is probably a direct consequence 
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of the mutual influence of dark and luminous components over the dissipative collapse of these 
objects. 
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